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ABSTRACT 

We develop an idealized inside-out formation model for disk galaxies to 
include a realistic mix of galaxy types and luminosities that provides a fair match 
to the traditional observables. The predictions of our infall models are compared 
against identical models with no-size evolution by generating fully realistic 
simulations of the HDF, from which we recover the angular size distributions. 
We find that our infall models produce nearly identical angular size distributions 
to those of our no-size evolution models in the case of a = geometry but 
produce slightly smaller sizes in the case of a f2 = 1 geometry, a difference we 
associate with the fact that there is a different amount of cosmic time in our 
two models for evolving to relatively low redshifts {z ^ 1 — 2). Our infall models 
also predict a slightly smaller (11% - 29%) number of large (disk scale lengths 
> 4:h'^Q kpc) galaxies at z ~ 0.7 for the CFRS as well as different increases in 
the central surface brightness of the disks for early-type spirals, the infall model 
predicting an increase by 1.2 magnitudes out to 2; ~ 2 (f2 = 0), 1 [Q = 1), while 
our no-size evolution models predict an increase of only 0.5 magnitude. This 
result suggests that infall models could be important for explaining the 1.2-1.6 
magnitude increase in surface brightness reported by Schade et al. (1995, 1996a, 
1996b). 



Subject headings: Galaxy Evolution - Galaxy Scale-lengths 
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1. Introduction 

Infall formation of disk galaxies predicts an inside-out formation process that produces 
galaxy sizes smaller than expected in models without intrinsic size evolution (Cayon, Silk 
& Chariot 1996). Prantzos & Aubert (1995) present a detailed comparison of several star 
formation models in an effort to explain the observed properties of the Galactic disk. The 
favoured model is one in which the star formation rate (SFR) is proportional to the gas 
surface density and to the differential rotation rate and in which gas evolution is driven 
by infall of unenriched gas during a certain period of time in the life of the galaxy (Wang 
& Silk 1994). In a first attempt to test infall formation of field disk galaxies, Cayon et 
al. (1996) noted that the observed trend toward small sizes at faint magnitudes and high 
redshifts is clearly reproduced by these models. 

In this Letter, we extend the development of the idealized infall formation model for 
disk galaxies by incorporating a range of different galaxy types and luminosities given by 
a slightly refined version of the Pozzetti, Bruzual, & Zamorani (1996) pure luminosity 
evolution (PLE) model. We aim to determine the extent to which these models differ from 
models with no-size evolution in the prediction of the angular size distributions, the number 
of large-intermediate redshift disks, and the evolution in surface brightnesses observed in 
disk galaxies at higher redshifts (Schade et al. 1995, 1996a, 1996b; Simard & Prichet 
1997). All calculations are performed with Ho — 50 km/sec/Mpc and expressed in Vega 
magnitudes unless otherwise noted. 



2. Description of Model 

Infall formation models have been favored by studies of the disk of our own Galaxy 
(Ferrini et al. 1994; Dopita & Ryder 1994; see Prantzos & Aubert 1995 for a detailed 
comparative analysis) and nearby galaxies (Ryder & Dopita 1994). A model for disk 
galaxies based on infall formation was presented in Cayon, Silk & Chariot (1996). Only a 
brief description is included here. A phenomenological model satisfying the observational 
requirements is one with infall in which the SFR depends on radius and time following a 
Schmidt-type law 

SFR{r, t)^{l- R)-'^^ Mopc-^Gyr-i, (1) 

where R 0.32 is the returned fraction of mass that was formed into stars, T.g{r,t) is the 
gas surface density at radius r and age t, and Tg{r) — [0.3(1 — R){rQ/r)]~^ Gyr is the gas 
consumption timescale {tq — 8.5 kpc). Normalizing the infall rate density of metal- free gas 
to the total surface density Ejot of stars plus gas observed at age T Gyr (corresponding to 
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z ~ 0) and using the instantaneous recycling approximation (IRA) , the evolution of the gas 
surface density is given by 

where t/ is the infall time scale. As in our previous work, for galaxies with Milky-way 
luminosities, i.e., M^,. fa —20, we take T,tot to equal the sum of the observed gas and the 
stellar surface density necessary to yield the assumed {z ^ 0) surface brightness profile. For 
galaxies of different luminosities, however, we take the star formation profile to be equal to 
the star formation profile for Milky- Way galaxies, but stretched by an appropriate factor in 
radius to produce the luminosity in question. 

Closely following the Pozzetti et al. (1996) PLE model, we develop a PLE model in 
terms of the above formalism for both an open and a flat universe, where we use three 
different morphological types: ellipticals, early-type spirals, and late-type spirals. To 
improve the agreement of our model with the number counts observed in the HDF for 
an open geometry, we increase the normalization and steepen the slope of our luminosity 
function slightly over the values given by Pozzetti et al. (1996). With these changes, we 
find that our models fit the number counts, redshift distributions, and luminosity functions 
as well as does the Pozzetti et al. (1996) model. 

We tune both the age T and the gas infall time scale Tf to match the z = colours 
and to provide a rough match with the star formation histories given for the Sabc and Sdm 
morphological types specified in the Pozzetti et al. (1996) PLE model with exponentially 
decreasing star formation. We take late-type spirals to have pure exponential profiles, and 
early-type spirals to have an exponential bulge added to this profile with a scale length 
that is 0.086 that of the disk (Courteau, de Jong, & Broeils 1996) with bulge-to-total ratios 
similar to those given in King & Ellis (1985). At redshift zero, we assume that our disks 
follow Freeman's law, i.e. = 21.65 mag/arcsec^ (Freeman 1970). For galaxies with 
Mhj > —18, however, we take the entire galactic population to be composed of late-type 
exponential disks with fi^j = 22.5 to account for the observed lack of early- type galaxies at 
low luminosities (Binggeli, Sandage, & Tammann 1988) and their lower surface brightnesses 
(McGaugh & de Blok 1997). For ellipticals, we adopt the size- luminosity relationship given 
by Binggeli, Sandage, & Tarenghi (1984). To perform the calculations for the described 
model, we break up our disk galaxies into approximately 30 annuli, on which we perform the 
spectral synthesis for each using a recent version of the Bruzual & Chariot tables (Leitherer 
et al. 1997). We summarize our choice of parameters in Table 1. 

We compare this model, hereinafter referred to as the "Infall" model, with a nearly 
identical "No Infall" model, in which there is no evolution in the scale length on which star 
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formation occurs, even though the star formation history for each morphological type in 
this model is the same as for the "Infall" model. Furthermore, we ensure that the surface 
brightnesses in the bj band for both models are exactly the same at 2; = 0. 

3. Evolution in Central Surface Brightness and Size 

For our first set of comparisons, we look at the extent to which the "Infall" models 
differ from "No Infall" models in terms of the evolution of different morphological types 
in the fio - tm plane. Because of the presence of the bulge in the center of the early-type 
spirals, we take the central surface brightness [Xo to equal the surface brightness extrapolated 
assuming an exponential disk profile from the surface brightnesses at ^ 0.5 and 2.0 
disk scale lengths-a procedure that is intended to be an approximate estimate of the 
central surface brightnesses that would be determined through the multicomponent fits 
(Schade et al. 1996a). Values in this plane are presented in Figure 1 for the models 
considered, from redshifts corresponding to an age of 1 Gyr down to 2; ~ 0. The values 
for L* Sab-Sbc and Scd-Sdm galaxies for the "Infall" model are shown as solid lines in 
panels (a) and (b), together with the predictions one magnitude above and below (shaded 
areas). Brightening of approximately 1.2 magnitudes occurs for early types, from z = to 
z ~ 2{VL = 0.0), 1{Q = 1.0), the rate of brightening being slower in the more open geometry 
due to the longer time between the formation of the galaxy and 2; ~ 1. These values would 
be in agreement with recent observations (Schade ct al. 1995, 1996a, 1996b) if early-types 
are dominating the samples. The faster the evolution in the central parts of the disk (larger 
T/t/ corresponding to earlier type disks) the larger is the amount of brightening observed. 
By contrast, our "No Infall" models, presented in panels (c) and (d) of Figure 1, show no 
such differential evolution, so consequently the early-type galaxies in these models produce 
only ~ 0.5 magnitude compared to the ~ 1.2 magnitude of brightening observed in the 
"Infall" models. These results suggest that inside-out formation scenarios, such as those 
we sketch in this paper, would not only make a difference in the observations relative to 
scenarios with no-size evolution, but they also could play an important role in the apparent 
surface brightness evolution observed to z ~ 1. 

As is clear from Figure 1, galaxies in our "Infall" models evolve not only in surface 
brightness but in size relative to "No Infall" models. In light of the present compilation 
of size information for the CFRS (Lilly et al. 1997), we make a simple estimate for the 
expected decrease in number of large galaxies predicted in the CFRS at intermediate 
redshift on the basis of our "Infall" models relative to our "No Infall" models. We find that, 
for the CFRS selection criteria (17.5 < Iab < 22.5), there are 807 (852) large (disk scale 
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lengths > 4 kpc) galaxies per redsliift interval per degree^ at z ~ 0.7 for our f2 = (f2 = 1) 
infall models compared to 719 (609) for our i7 = (i7 = 1) "No Infall" models, a decrease 

of 11% {n = 0) to 29% {n = 1). 

Finally, we examine the differences in the angular size distributions produced by the 
models in the context of the HDF, the deepest and highest resolution optical observation 
to date. To this end, we generate detailed images with areas equal to roughly 3.8 times 
the area of the three WF chips and with properties which closely match those observed in 
the F814W HDF images. In our simulations, we place galaxies with random orientation 
angles and no extintion on the mock images assuming no correlation, add Poissonian noise, 
smooth over a Gaussian kernel with a 0.06-arcsec sigma, and then degrade the image by 
adding some noise. After adding this noise, we smooth the images again over a Gaussian 
kernel with a 0.02-arcsec sigma to reproduce approximately the drizzled noise properties of 
the images (Williams et al. 1996). We use SExtractor 1.2b5 (Bertin & Arnouts 1996) to 
make a catalogue and to determine both the magnitude (SExtractor's MAG_BEST) and the 
half-light radiiQ of the objects in the F814W images for each of our models as well as the 
F814W images (version 2) of the three WF chips (1634 objects total). For a detection, we 
require an object, after smoothing with the PSF, to be composed of at least 10 contiguous 
pixels at 2a above the background. 

As was our goal, in Figure 2 we demonstrate that the f2 = models provide fair fits 
to the HDF number counts, decomposed into three different bins of half-light radii (though 
we are 20-40% low in the faintest magnitudes bins). We then compare the angular size 
distributions recovered from the simulations for our "Infall" models with those recovered 
from both our "No Infall" models and the HDF in Figure 3, where we decompose our 
angular size distribution into its constituent morphological types. We find that while both 
our n = models are in fair agreement with the angular size distributions found in the 
HDF, they also appear to be in fair agreement with each other, even for the early-type 
galaxies which evolve so very differently in the two models (Figure 1). For our f2 = 1 
models, however, there is a clear difference between the angular size distributions for the 
"Infall" and "No Infall" models, particularly among the early-type spirals. 

To understand this difference, it is important to note that the profiles of galaxies in 
our "Infall" models tend to approach those in our "No Infall" models over time, and so 
the convergence begins in the center and grows outward. Fundamentally, both our = 
models produce similar distributions of recovered half-light radii because at any redshift 



^We take the half-light radius to equal the radius of the aperture which contains half the light as 
determined by SExtractor's best estimate of its total (MAG_BEST). 
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only those central parts of the galaxy for which the profiles in the two models have started 
to converge can be seen, these radii being small at high redshifts and large at low redshifts 
as a result of (1 + z)'^ cosmological dimming. Galaxies in our Q = 1 models, however, have 
less time to evolve to relatively low redshifts where cosmological dimming is less significant 
(note that in Figure 1 the redshift corresponding to the peak in surface brightness is lower 
for our Q = 1 models); therefore, the apparent differences between the profiles produced 
by the two models can be more readily observed. In Figure 3, we plot the distribution of 
true half-light radii (calculated from our code which placed each galaxy on our simulated 
images) for the two Q = models to illustrate the reality of the actual but relatively 
unobservable differences between these two models. 

In summary, we have extended the idealized infall model for disk formation to include 
a realistic mix of galaxy types and luminosities which provide a rough match to the number 
counts, redshift distributions, and angular size distributions in the HDF. We compare the 
predictions of these models against identical models with no-size evolution to see how our 
infall models would affect the observations. For these infall models, we find that the number 
of large disks (disk scale lengths > 4 kpc) is ll%-29% smaller than for our no-size evolution 
models at 2; ^ 0.7 for the CFRS. We also find increases in the central surface brightness of 
disk galaxies of about 1.2 mag/arcsec^ out to 2; ~ 1(^2 = 1), 2(^2 = 0) for early-type spirals, 
very different from the 0.5 mag/arcsec^ increase found for models with no-size evolution. 
Our finding strongly suggests that inside-out formation scenarios for disk galaxies should 
have observable consequences at high z and could be important for understanding the 
observed increase in surface brightness of disks out to z 1. 

For our Q = models, we find no significant differences between the recoverable angular 
size distributions predicted on the basis of our infall models and those predicted on the 
basis of our models with no-size evolution. We suspect that significant differences between 
these models would have been obtained if the early-type galaxies had been allowed to begin 
forming at lower redshifts, as we found in the case of our f2 = 1 models. However, in this 
case, there would have not been as many galaxies at high redshifts and our models would 
have undcrpredicted the number of galaxies at faint magnitudes. We therefore suspect that 
for any simple single-formation-epoch PLE models which fit the number counts, infall and 
no-size evolution models will tend to produce very similar angular size distributions given 
the current state of observations. Therefore, if our infall models are correct, the fact that 
other workers (Ferguson & Babul 1997) have found angular size distributions with no-size 
evolution PLE models which tend to agree with the HDF observations is not very surprising 
in fight of the present result. 

We would hke acknowledge Ian Small for his suggested procedure for the determination 
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Fig. 1. — Disk central surface brightness /ig versus half-light radius rhi for Sab-Sbc (top lines) 
and Scd-Sdm (bottom lines) galaxies. The redshift interval considered in the calculations 
extends from z = up to a redshift that corresponds to the time at which galaxies have 
been evolving 1 Gyr in the different models. Solid lines denote the values corresponding to a 
M* galaxy with the shaded areas covering one magnitude above and below this. Predictions 
of the "Infall" models are presented in panels (a) and (b) while the "No Infall" models are 
presented in panels (c) and (d). Predicted values in = 0.0 and Q = 1.0 cosmologies appear 
in the left panels and right panels respectively. Tick marks denote the locations at different 
redshifts in steps of 0.2 from z — to z — 1 and in steps of 1 from z — 1 to the highest 
redshift indicated for each model. 
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Fig. 2. — Comparison of the number counts recovered from the HDF with those recovered 
from the simulations. The observed number counts are broken up in terms of their measured 
half-hght radii into three bins Vhi < 0.18 arcsec, 0.18 arcsec < rhi < 0.30 arcsec, and 0.30 
arcsec < rhi- The error bars are the one-sigma poissonian values. 
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Fig. 3. — Comparison of the half-light radii distributions recovered from the HDF 
(histogram) with one-sigma error bars against those recovered (r/j/(rec)) from our simulations 
in the magnitude bins (23 < mpsuwiy^S^) < 25; 25 < mpsuwiy^S^) < 27). The angular 
size distributions recovered from the simulations are decomposed into morphological types. 
The true half-light radii (jhi^in)) distributions of higher rcdshift {z > 1) objects recovered 
in the magnitude bin 25 < m^8i4VK(Vega) < 27 are compared to the recovered distributions 
to illustrate both the errors in our procedure for determining the half-light radii and the 
fact that differences in half-hght radii between the "Infall" and "No Infall" models are more 
apparent in the true half-light radii distributions than in the recovered distributions. 
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Table 1. Model parameters" 



Q 


Type 


a 


^6. 


0S 








SFR 


T(Gyr) 


0.0(1.0) 


E/SO 


-0.48 


-20.87 


0.95 


1 








16.0 (12.7) 




Sab-Sbc" 


-1.32 


-21.14 


0.69 


0.2 


21.65 




7.2 Gyr(3.9 Gyr) 


18.0(12.7) 




Sab-Sbc" 


-1.32 


-21.14 


0.69 


0.13 


21.65 


^/ = 


7.2 Gyr(3.9 Gyr) 


18.0(12.7) 




Scd-Sdm" 


-1.32 


-21.14 


0.76 





21.65 


= 


- oo Gyr(oo Gyr) 


18.0(12.7) 




Scd-Sdm-^ 


-1.5 


-21.14 


1.35 





22.5 




: OO Gyr (oo Gyr) 


18.0(12.7) 



^Scalo IMF assumed for the E/SO and Sab-Sbc types and Salpeter IMF for Scd-Sdm types. 
'^Bulge-to-total luminosity ratio in hj band at 2; = 0. 
■^Exponential SFR characterized by decay time t\ = 1 Gyr. 
'^l Gyr burst SFR. 

^LF truncated to include only galaxies with M^j < —18. 
^LF truncated to include only galaxies with M^j > —18. 



